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Abstract 

A series of active flow control experiments were 
recently conducted at high Reynolds numbers on a 
generic separated configuration. The model simulates 
the upper surface of a 20% thick Glauert-Goldschmied 
type airfoil at zero angle of attack. The flow is fully 
turbulent since the tunnel sidewall boundary layer 
flows over the model. The main motivation for the 
experiments is to generate a comprehensive data base 
for validation of unsteady numerical simulation as a 
first step in the development of a CFD design tool, 
without which it would not be possible to effectively 
utilize the great potential of unsteady flow control. 
This paper focuses on the dynamics of several key 
features of the baseline as well as the controlled flow. 

It was found that the thickness of the upstream 
boundary layer has a negligible effect on the flow 
dynamics. It is speculated that separation is caused 
mainly by the highly convex surface while viscous 
effects are less important. The two-dimensional 
separated flow contains unsteady waves centered on a 
reduced frequency of 0.8, while in the three 
dimensional separated flow, frequencies around a 
reduced frequency of 0.3 and 1 are active. Several 
scenarios of resonant wave interaction take place at the 
separated shear-layer and in the pressure recovery 
region. The unstable reduced frequency bands for 
periodic excitation are centered on 1.5 and 5, but these 
reduced frequencies are based on the length of the 
baseline bubble that shortens due to the excitation. 
The conventional swept wing-scaling works well for 
the coherent wave features. 

Reproduction of these dynamic effects by a 
numerical simulation would provide benchmark 
validation. 

Nomenclature 

steady blowing momentum coefficient, = J j cq 

<c M > oscillatory blowing momentum 

coefficient, = < 7 ? > fcq 

Cy combined blowing momentum 

coefficient, = (c p ;< c p >) 
c model chord 

C p wall pressure coefficient, = (P - P s )/ q 


f oscillation frequency [Hz] 

F + reduced frequency, = / * X fU ^ 

F fundamental excitation frequency 

h slot height or width 

H height of separated region 

2 

J momentum at slot exit, = phU j 

M Mach number 

nF integer multiple of the excitation freq 

P pressure 

q free stream dynamic pressure, q - 1 12 pU^ 
R reattachment 

R c chord Reynolds number, R^U^civ 

S separation 

T temperature 

U, u average and fluctuating streamwise velocity 
U<p phase velocity 

x/c normalized streamwise location 

X sp distance from actuator or separation to 

reattachment 

z spanwise location 

A sweep angle [deg] 

v kinematic viscosity 

p density 

0 BL momentum thickness 

Abbreviations 

BL boundary layer 

LE leading edge 

TE trailing edge 

< > phase locked values 

Subscripts 

b baseline parameters 

c cavity 

f related to excitation frequency 

j conditions at blowing slot 

e conditions at BL edge 

oo free-stream conditions 

2D two-dimensional 

3D three-dimensional 

Superscripts 

' root mean square of fluctuating value 
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1. Introduction 

This paper is the third in a series of publications 
describing experiments that were aimed at improving 
our understanding of controlling separated flows at 
flight Reynolds numbers and providing a 
comprehensive database for validation of unsteady 
CFD design tools. Specifically, effects of the upstream 
boundary layer thickness, compressibility, sweep and 
location for introduction of the control input are 
explored. A previous publication 1 presented some of 
the experimental results for incompressible 2D flow 
over the same model. Periodic excitation is the primary 
control mechanism, but steady mass transfer is also 
used, for comparison of effectiveness and also in order 
to provide fully attached steady flow, since the baseline 
uncontrolled flow contains a large turbulent separation 
bubble. 

The model simulates the upper surface of a 20% 
thick Glauert-Goldschmied type airfoil and has some 
features in common with a backwards-facing step with 
the exception that the flow on the model can be fully 
re-attached with effective control. Without control, the 
flow separates at the highly convex region of the model 
(x/c=0.65, Fig. 1) and a large turbulent separation 
bubble is formed. Periodic excitation is applied from 
the slot located at x/c=0.64 or at x/c=0.59 to gradually 
eliminate the separation bubble. The effects of steady 
suction or blowing were also tested. The dynamics of 
these steady flow control methods are compared to 
those of the periodic excitation. 

A brief summary of the incompressible 2D 
findings that were presented in Ref. 1 is repeated here 
for completeness. The baseline flow was carefully 
documented and found to be weakly dependent on 
Reynolds number. The baseline flow is fully turbulent 
so laminar-turbulent transition, due to the active 
separation control, does not effect the data trends. The 
thickness of the upstream boundary layer was 
monitored and controlled. A reduction of about 40% in 
the boundary layer momentum thickness had only a 
minor effect on the baseline separation and its control. 
The spanwise uniformity of the wall pressure 
distribution was found to be very good and improved 
as the separation was controlled. The baseline 
reattachment point was found to be about one 
separation height downstream of the peak in the wall 
pressure fluctuations. The separation height is taken as 
the model elevation at separation (i.e. x/c=0.65, 
H=0.115c). The level of this peak increased as 
reattachment moved forward due to control, regardless 
of the control method. Steady suction or blowing with 
a momentum coefficient of 2-4% is required to fully 
reattach the flow to the model surface and recover the 
ideal pressure distribution. Active control using 
periodic excitation is comparable to steady suction and 
significantly more effective than steady blowing. 
Periodic excitation is capable of reattaching the flow 


but is not capable of reproducing the same pressure 
jump across a slot as suction does. This is presumably 
since the periodic excitation relies on enhanced mixing 
that is by nature a convective phenomena rather than a 
local one that is generated by the severe suction. It was 
found that the superposition of weak suction on the 
oscillatory excitation significantly enhances the 
receptivity of the separated shear layer to the 
fundamental excitation frequency. Detailed dynamics 
of the wall pressure fluctuations that explain some of 
the observed features will be discussed. 

The measurement tools used for the dynamic 
testing are 12 dynamic pressure transducers that are 
positioned next to the static pressure taps (x symbols 
on Fig. 2). The dynamic pressure measurement allows 
the evaluation of random and coherent pressure waves. 
In addition, the measurement of the cavity pressure 
fluctuations provides an opportunity for a 
comprehensive description of the control input. 

The dynamics of pressure fluctuations measured 
over a wide range of separated flow configurations 
were reviewed in Reference 2. It was shown that the 
level of pressure fluctuations increase downstream of 
separation, and peaks close to and upstream of 
reattachment, at values of 4-10% the dynamic pressure. 
The dimensionless dominant frequency, based on the 
length of the separated region and on the free stream 
velocity was between 0.5 and 0.8, The excellent review 
of impinging shear-layers by Rockwell 2 3 , supports the 
above findings and discusses linear as well as non- 
linear evolution of the separated and reattached shear 
layer. 

Two extensive direct numerical simulations (DNS) 
of separated flow were recently published 0 . For a 
backwards-facing step, the separated free shear-layer 
was unsteady with a dominant range of F* of one at 
separation and it reduced to 0.5 at reattachment. The 
reduction in F* was an indication of dissipation of the 
small structures or vortex merging. The F* was based 
on the length of the bubble, X sp ~6.3 step heights. The 
phase speed of the F*=0.5 structures increased 
downstream of reattachment. The calculated wall 
pressure fluctuations peaked just upstream of 
reattachment, with a value lower than experimentally 
measured. The second DNS 5 deals with separated and 
re-attached flow due to an adverse pressure gradient 
imposed on a flat plate boundary layer. The main 
difference between this flow and the backwards-facing 
step is the lack of a fixed separation point, that indeed 
was found to fluctuate significantly in the streamwise 
direction. The dominant F^ as well as the peak of C p ’ at 
reattachment were similar in both flows. 

This paper does not pretend to provide a 
comprehensive description of the flow field dynamics, 
as this could not be provided using only the wall 
pressures. Rather it attempts to deepen the 
understanding of several key features and provide 
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validation cases for future CFD. A brief review of the 
experimental set-up is provided in Sec. 2. The 2D 
baseline flow is discussed in Sec. 3. The effect of the 
upstream BL thickness is discussed in Sec. 3.2. The 
effect of sweep on the baseline flow spectra is 
discussed in Sec. 3.3. The remainder of the paper deals 
with the controlled flow. Sec. 4.1 describes the effects 
of the excitation frequency and its amplitude. The 
effects of sweep on the 3D controlled flow are 
discussed in Sec. 4.2. Finally, the unsteady effects of 
weak steady suction and periodic excitation are 
compared in Sec. 4.3. 

2, Description of the Experiment 

2.1 Overview 

This section contains a brief description of the 
experiment. A more detailed description can be found 
in References 1 and 6. 

2.2 Model 

The 200 mm chord model simulates the upper surface 
of a 20% thick airfoil that is a variation on the Glauert 
Glas II airfoil 7 . A moderate favorable pressure gradient 
up to 55% of the chord is followed by a severe adverse 
pressure gradient, imposed by the highly convex 
surface at x/c-0.6, that relaxes towards the trailing 
edge (see Fig. 1). The model was mounted on the 
sidewall of the tunnel and therefore the flow over the 
model was fully turbulent. Two alternative excitation 
slot locations, x/c=0.59 and x/c=0.64, were available. 
The slots were about 0.25% chord wide (0.50 mm 
±20%), and allowed an almost tangential downstream 
introduction of momentum (the slots are inclined at 30° 
to the surface because of manufacturing considerations, 
see Fig. 1). The model design also enabled testing at 
two sweep angles, A=0 degrees and A=30 degrees. A 
top view of the 2D configuration is shown in Fig. 2. 
The symbols in Fig. 2 show the distribution of the 
static pressure taps and dynamic pressure transducers, 
while the vertical thick lines mark the location of the 
end plates that isolate the boundary layers on the floor 
and ceiling of the wind tunnel from affecting the flow 
over the model. Ten dynamic pressure transducers 
(marked by x symbols on Fig. 2) are installed under the 
model surface, inside small volume cavities. The 
cavities are connected to the surface of the model by 
tiny orifices, 0.254 mm in diameter. The effect of this 
installation on the frequency response of the dynamic 
pressure transducers was studied using a 
comprehensive bench-top calibration and occasional 
in-situ testing by a Piezo-electric actuator. In addition, 
a transducer was flush mounted near the trailing edge 
of the turntable (x/c=1.4) next to a recessed pressure 
transducer and their readings are compared. The full 
scale of the dynamic pressure transducers is 10 psid. 
They are referenced to the static wall pressures 
immediately next to the orifice locations or to the wind 
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tunnel plenum pressure in order to maintain optimal 
resolution even at static pressures that exceed the 
transducer’s range. One dynamic pressure transducer is 
also installed inside the model cavity, midway between 
the end plated and about 30-mm from the slot exit. It is 
used to monitor the cavity pressure oscillations and to 
correlate the wind tunnel experiments with the bench- 
top tests that were conducted in order to correlate the 
slot exit velocities with the cavity pressure 
oscillations 1 . 

2.3 Oscillatory Blowing System 

The oscillatory blowing system is capable of 
introducing a wide range of steady and periodic 
momentum combinations from the cavity inside the 
model to the external flow. More details can be found 
in References 1 and 8. 

2.4 Boundary Laver Measurement System 

The boundary layer probe consisted of a Pitot-probe 
that was connected to both a static and a dynamic 
pressure transducer. The probe was mounted to the 
backside of the model inside a temperature-controlled 
box. The boundary layer probe penetrated to the flow 
side of the model through six alternative slots in the 
turntable. The, internal diameter of the probe tip was 
0.254mm and when touching the wall, the probe 
integrated the total pressure from 0.1mm to 0.35mm 
from the wall. A more detailed description of the 
boundary layer measurement system can be found in 
Reference 1. 

2.5 Experimental Uncertainty 

The experimental uncertainty in the determination of 
C p 7 amplitude was estimated to be ±15%. The 
uncertainty of the phase of the fundamental excitation 
frequency and its second harmonic was estimated to be 
±10 deg. This leads to an uncertainty of the calculated 
phase velocities of ±6%. The estimation of the <c M > 
was within ±25% of the quoted values, but this value is 
probably lower for most of the non-cryogenic 
conditions cited in this paper. The steady mass flux 
momentum coefficients are within ±5% of the cited 
values or 0.01% absolute value (the bigger of the two). 
Additional uncertainty information can be found in 
Refs. 8 and 1 . 

2.6 Experimental Flow Conditions 

The experiments were conducted at Mach numbers 
from 0.20 to 0.70 and chord Reynolds numbers ranging 
from 2.4x1 0 6 to 39x1 0 6 . 
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3. Dynamics of the Baseline Flow 

3.1 Baseline 2D Flow 

The Reynolds number has a weak effect on the model 
pressure distributions, as seen from the 2D data in Fig. 
3. Note the C p ’ ordinate is on the right side of the 
figure. The flow separation (at x/c=0.65) is not affected 
by the presence of the blowing slot, by an interaction 
between the external flow and the cavity 1 , or even by 
sweeping the model to A=30 deg (Fig. 3 here and Fig. 
3 of Ref. 6). Mean and fluctuating pressure 
distributions over a R c range of 7 to 26x1 0 6 were 
compared (Fig. 5, Ref. 1), showing very weak R c 
sensitivity. Fig. 3 compares data measured at 
R c =2.4xl0 6 with data taken at R c =16xl0 6 . The very 
small effect the Reynolds number has on the wall 
pressure distributions facilitates studying the pressure 
field dynamics at relatively low R c (i.e., non-cryogenic 
temperatures), where the flat range of the installed 
dynamic pressure transducers’ frequency response 
extends up to 2.5 kHz. 

The fluctuating wall pressures (C p r ), also shown in 
Fig. 3, are very low upstream of separation. A 
significant increase in C p is seen for x/c>0.64 (the 
thick dashed vertical line indicates the slot location) 
and a further increase as reattachment is approached 
(x/c=l.l). Reattachment is assumed to be between the 
location of C p ' max (denoted by x/c(C p f maA )) and one 
separation height downstream of x/c(C p ' max ) 2 ’ 4,5 . The 
separation height is taken as the model elevation at 
separation (i.e. x/c=0.65, H=0. 1 15c). The wall pressure 
fluctuations decrease downstream of reattachment as 
the turbulent boundary layer readjusts to the zero 
pressure gradient flow. The pressure amplitude spectra 
at R c =4.2xl0 6 are presented in Fig. 4. Note that the 
spectra of x/c>0.85 are shifted up by a decade with 
respect to previous x/c spectra. Very close to 
separation (x/c=0.73), a broad band of active 
frequencies centered on F^O.S can be observed to 
slightly rise above the background random pressure 
fluctuations. Note that the frequencies were reduced 
using the length of the separated region (~c/2) and the 
free stream velocity. These frequencies are probably 
related to a Kelvin-Helmholz type of instability that is 
independent of the Reynolds number or to the laminar- 
turbulent state of the boundary layer prior to 
separation. The low frequency range fills very quickly 
above the bubble (x/c=0.85), and the overall pressure 
fluctuation level increases. A further increase in C p \ 
accompanied by a halving of the dominant frequencies 
to F*=0.4, occurs at reattachment (x/c=1.11). The low 
frequencies decay very rapidly in the attached 
boundary layer (x/c=1.42), leaving a dominant band of 
F* centered on F*=0.5. It is highly probable that sub- 
harmonic resonance takes place above the bubble, but 
the resolution of the measurements did not allow a 
clear identification of the process dynamics. This point 
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will be addressed again in Section 4.1 discussing the 
controlled flow. 

3.2 The Effect of the Upstream Boundary Laver 
Thickness 

It was demonstrated 1 that the thickness of the upstream 
boundary layer has very little effect on the baseline and 
controlled pressure distribution in the range R^-7- 
26xI0 6 , Furthermore, it was demonstrated that thinning 
of the upstream boundary layer by distributed suction 
upstream of the model by as much as 43% (based on 
the momentum thickness) does not have a measurable 
effect on the integral parameters. This point is 
addressed again here at a lower chord Reynolds 
number (R c =4.2xl0 6 , M=0.25, 2D flow, and x/c=0.64 
slot) where the effect is expected to be larger than at 
high Rc. Fig. 5 shows the mean total pressure profiles 
measured at 1.25c upstream of the model leading edge. 
The strong suction reduced the momentum thickness 
by 35% and the shape factor from 1.29 to 1.17. The 
Reynolds numbers based on the upstream BL 
momentum thickness (Re©) were 3 1 x 1 0 3 and 20xl0 3 
for the baseline and boundary layer control (BLC) 
flows respectively. The corresponding model pressure 
distributions are presented in Fig. 6. The main effect of 
the thinner upstream boundary layer is a stronger 
deceleration of the flow upstream of the LE, in 
accordance with the ideal flow prediction. The pressure 
distributions are identical from x/c=0 to separation, 
x/c=0.65. The thinner boundary layer flow separates at 
a slightly higher C p and generates a lower C p ’ at 
x/c=0.67 (right hand ordinate). The pressure recovery 
downstream of the bubble is similar, but the 
reattachment C p is higher for the thinner boundary 
layer flow. This is accompanied by a higher C p \ It is 
speculated that the higher C p ’ at reattachment causes 
the higher C p that is measured downstream of 
reattachment, due to the impingement of the large 
coherent structures on the wall. The pressure amplitude 
spectra for x/c=0.73 are presented in Fig. 7. The first 
thing to note is the very weak effect of the thinner 
boundary layer on the spectra. The thinner boundary 
layer flow is characterized by a reduced level of 
turbulent activity in the range 0.3<F + <2. At 
reattachment (data not shown), x/c= 1 . 1 1 , the thinner 
boundary layer flow exhibits higher C p ’ at F*>1, while 
downstream of reattachment, at x/c=1.42, the thinner 
boundary layer exhibits slightly lower C p ’ at F^cO.S. 
Again, the overall very small effect the thinner 
boundary layer has on the spectra of the wall C p ’ is the 
main finding here. 

3.3 The Effect of Sweep on the Baseline Spectra 

It was shown 6 that sweep has two major effects on the 
model wall pressures: first, the bubble is shorter and 
second, the level of C p ’ from separation to 
reattachment is higher. Certainly, these two features are 
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interconnected. The higher level of wall pressure 
fluctuations is a manifestation of increased turbulence 
activity, in the form of large coherent structures that 
mix more effectively across the separated shear-layer, 
causing earlier reattachment. Fig. 3 presents a 
comparison between the 2D and 3D wall pressures at 
M=0.25 and R c =2.4xl0 6 and 4.2x1 0 6 . These data are 
essentially identical to the high Rc data shown in Fig. 9 
of Ref 6. The power spectra for A=0 and A=30 are 
compared in Figs. 8a through 8d. The overall level of 
C p ’ immediately downstream of separation (x/c=0.73, 
Fig. 8a) increases significantly due to sweep. Instead of 
the single, wide band of C p ‘ centered on F^O.8 in the 
2D flow, two peaks can be seen in the 3D flow. A wide 
and flat ridge, centered on F + =l, and a very distinct 
peak centered on F + =0.3 are evident. It is assumed that 
the additional peak at F*=0.3 is due to sweep, and 
might even be set-up dependent, due to the finite span 
of the model. The F^l peak scales favorably with the 
reduced length of the bubble (X sp was maintained as 
c/2, while it actually is shorter than that of the 2D 
flow). Above the bubble (i.e. x/c=0.85, Fig. 8b), the 
F 4 ^! peak shifts to lower frequencies while the peak at 
F^O.3 is still present. The spectra of the 2D flow at 
x/c=0.85 are almost featureless. At reattachment 
(x/c= 1 . 1 1, Fig. 8c), the F* =0.3 peak is still seen for the 
3D flow while the overall level of C p ' decreases. The 
2D flow shows a sub-harmonic peak at F* =0.4, while 
there is no evidence of sub-harmonic resonance in the 
3D flow. Downstream of reattachment (x/c=1.42, Fig. 
8d) the C p ’ spectra for both the 2D and 3D flows 
decrease over the entire F* range, but the decrease for 
the 3D flow is significantly faster, since x/c=1.42 is 
further downstream from reattachment for the 3D flow. 

4.0 Dynamics of the Controlled Flow 

4.1 The Effect of Excitation Frequency in 2D Flow 
Fig. 9 presents 2D pressure distributions measured at 
R<.=4.2xl0 6 and M=0 .25 with control applied from the 
x/c=0.64 slot at F* =0.5 and F + =l. Even though <c M > at 
F + =l is only a third of <c p > at F^O.5, the mean 
controlled pressures upstream of the slot are identical. 
A small bubble remains for the controlled flow using 
F + =1 , while flow separation is delayed using F^O.5 
with <c p >=0.3%. The total rms of the wall pressure 
fluctuations (C p ’) of the F*=0.5 and F*=l controlled 
flows at x/c=0.67 are identical. While C p ’ at F*=0.5 is 
amplified up to x/c=0.8, C p ’ at F + =l saturates at 
x/c=0.73 and begins to decay rapidly in the smooth 
pressure recovery region of the flow excited by F + =l. 
The mean pressures are identical for x/c>l.l. The 
approach currently used to determine <c M >, is to 
measure the cavity pressure fluctuations (Fig. 1) and 
correlate these with the measured velocity fluctuations 
in the absence of an external stream 1,8 . This procedure 
allows the evaluation of <c p > over the entire frequency 
domain of interest, not only at the excitation frequency 


and does not rely on the hot-wire de-rectification 
process that can distort the frequency content. Fig. 10a 
presents the <c p > for F^ =0.5 and F*=l that were 
calculated from the measured cavity pressure 
fluctuations. The fundamental excitation frequency 
contains 81% and 86% of the total excitation 
momentum for the fundamental of F + =0.5 and F*=I 
respectively. Note that <c p > of the 2 nd harmonic of 
F^O.5 is smaller by an order of magnitude than <c M > 
of the fundamental of F*=0. 5. The 2 nd harmonic of 
F*=l is smaller than the fundamental <c p > by almost 
three orders of magnitude. Comparing the total 
fluctuating momentum contained in the fundamental 
and its higher harmonics, in the range 0.5<F f <4, 
indicates that it is 0.38% and 0.14% for F*=0.5 and 
F*=I respectively, maintaining a ratio of about 3:1 
between <c p > of F*=0.5 and <c p > of F + =l. Fig. 1 0b 
compares the baseline and controlled pressure spectra 
at x/c=0.73, immediately downstream of the excitation 
slot and very close to the location of the baseline flow 
separation. The pressure amplitudes of both excitation 
frequencies are identical at this location, even though 
<c p > at F*=0.5 was three times higher at the excitation 
slot (Fig. 10a). The C p ’ of the 2 nd harmonic of F*=1 is 
similar to the C p ’ of the 3 rd harmonic of F f =0.5. These 
frequencies generate F^’s of 1.5 to 2, which were 
shown to be the most effective frequencies for 
separation control over this geometry at high Reynolds 
numbers 1 . It is demonstrated that the flow is most 
receptive to these F*’s, because the ratio between the 
excitation momentum and the measured C p ’ at these F* 
are the highest. Fig. 10b also shows that the overall 
level of the pressure fluctuation increases due to 
excitation, with the exception of F*<0.1, with respect 
to the baseline C p \ The spectra of the wall pressures in 
the recovery region (x/c=0.85) are presented in Fig. 
10c. The most striking feature is the strong decay of 
the higher harmonics of the excitation at F*>1, that 
were very active above the bubble (Fig. 10b). While 
F^O.5 maintains its momentum, F + =l decayed 
between x/c=0.73 and x/c=0.85, and the amplitude of 
the 2 nd harmonic of F^O.5 increased. Note that F*=1 
developed two side bands at F*~0.5 and at F^l.S, 
which is an indication of sub-harmonic resonance. The 
C p ’ downstream of reattachment are presented in 
Fig.lOd. The controlled flow is steadier than the 
baseline for F*<0.5, for both excitation frequencies. 
The decay rate of the F + =l related C p ’ is faster than the 
decay rate of F*=0.5 related C p \ The random pressure 
fluctuations of the flow that was excited by F + =0.5 are 
lower than the random C p ’ of the flow excited by F*=l . 
It will be shown later that this feature is related to the 
excitation level when using low F \ The C p ’ spectra at 
the end of the measurement domain are presented in 
Fig. lOe (x/c=1.42). For F*=l only the C p ’ at the 
excitation frequency and its 3 rd harmonic are higher 
than the baseline C p \ while for F* =0.5 the fundamental 
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and its two higher harmonics rise above the baseline 

cy. 

Additional information regarding the dynamics of 
the control process can he gained by studying the 
evolution of the excitation related pressure fluctuations 
(Cp’f) and the phase speed of the coherent pressure 
fluctuations. C p ’ f can be extracted from the C p ’ spectra 
for highly coherent flow or from the ensemble 
averaged pressure fluctuations, phase-locked to the 
excitation frequency. Fig. 1 la compares the C P V at the 
excitation frequency for the data presented in Figs. 9 
and 10. It is demonstrated that the amplitudes at 
x/c=0.67 are similar even though there was a factor of 
three in the input excitation momentum. The C p f at 
F^O.5 is amplified while at F*=l the C p ’ f saturates at 
x/c=0.73. The decay rate of F*=l is faster than the 
decay rate of F*=0.5. Similar trends were observed at 
high R c (Fig. 14c, Ref. 1). The phase velocities of the 
coherent wall C P ’ f were calculated using a 3 rd order 
polynomial fit to the evolution of the phase with x/c 
and normalized by the local velocity at the BL edge, 
U e , that was calculated from C p . The phase velocities 
above the bubble (Fig. lib, x/c~0.7) are about 0.3U e , 
close to the linearly expected values for low amplitude 
disturbances in attached boundary layers, and smaller 
than those expected in a 2D turbulent mixing layer in 
the absence of a pressure gradient. The phase velocities 
increase as reattachment is approached, and the phase 
velocity of F + =l accelerates faster. Downstream of 
reattachment both phase velocities approach 0.8, which 
is an acceptable value for turbulent boundary layer 
structures. 

The relationship between the fundamental 
excitation frequency and its second harmonic was 
evaluated for F^=0.5 . Fig. 12a presents the total rms of 
the wall pressure fluctuations (C p ’) as well as the level 
of Cp’, at F*=0.5 and at its second harmonic, (i.e. C p2 f , 
F^l). Fig. 12b presents the phase velocities of F*=0.5 
and its 2 nd harmonic. Note that the excitation input is 
contaminated by higher harmonics since C p ’f contains 
only 74% of the total C p ’ at x/c-0.67 (it contained 81% 
of <c p > at the excitation slot. Fig. 10a). The C p2 ( is 
about 13% of the total C p ’ at x/c=0.67, and it initially 
decays. The phase velocity of the fundamental (Fig. 
12b) indicates that it is a convective shear-wave, while 
the phase velocity of the 2 nd harmonic of the excitation 
frequency indicates that it is initially a standing wave 
(i.e. U^— »0). Between x/c=0.77 and x/c=1.25 the phase 
velocities of F and 2F are similar, within the 
experimental uncertainty. It is also noted that as F 
starts to decay, 2F amplifies (Fig. 12a) over the 
streamwise extent that corresponds to the pressure 
recovery region (i.e. 0.77<x/c<l, Fig. 9). These 
findings indicate that a resonant wave interaction is 
active in the reattachment region. In the attached 
region of the boundary layer (x/c=1.25), F continues to 


accelerate to 0.8U e , while 2F, whose amplitude decays 
slower than the amplitude of F also accelerates slower. 

4.2 Excitation Effect in 3D Flow 
It was demonstrated 6 * that the effectiveness of periodic 
excitation is not reduced at mild sweep angles. It was 
further shown that the integral parameters of the 
controlled 3D flow scale according to simple swept 
wing scaling 9 . Presently, the dynamics of the 3D 
controlled flow are examined and several key features 
of 2D and 3D controlled flows are compared. 

Fig. 13 presents baseline and controlled pressure 
distributions measured at R c =4.2xl0 6 , M=0.25 with the 
model swept to A=30 deg and an excitation slot 
located at x/c=0.64. The baseline reattachment is 
assumed at x/c~l.l and X sp =c/2 is used, for the sake of 
simplicity and ease of comparison with the 2D flow. 
The data indicate the effectiveness of F*=l with 
<c p >=0.04% due to the ability of this excitation mode 
to shorten the separation bubble. C p of F + =l decays 
immediately downstream of the excitation slot. The 
flow is less receptive to F^O.5, as indicated by C p of 
F*=0.5 and <c p >=0.04% at x/c=0.67 with respect to 
C p ’ of F*=l at the same <c M >. An increase of the 
excitation momentum of F^O.5, from <c p >=0.04% to 
<c p >=0.3%, results only in a minute increase of C p ’ at 
x/c=0.67. This could be partly due to the effective 
delay of separation, as indicated by the C p of the 
F^O.5, <c m >=0.3% excited flow. Separation delay, 
which is the purpose of the control strategy, alters the 
mean flow and its instability characteristics, making the 
attached flow less receptive to excitation, leading to 
saturation of the control effectiveness. 

The excitation momentum inputs, at F*=0.5 and 
F f =1.0 for low <c M > controlled flow data that are 
presented in Fig. 13, were calculated from the cavity 
pressure fluctuations and are plotted in Fig. 14a. The 
values of <c p > at the two fundamental excitation 
frequencies are identical. The higher harmonics of 
F*=0.5 (dashed line) are at worst three orders of 
magnitude smaller than the <c M > of the fundamental. 
The higher harmonics of F*=l are about two orders of 
magnitude smaller than the <c p > of the fundamental, 
with the 3 rd harmonic (i.e. F* ~3) larger than the 2 nd 
harmonic (F f =2). The wall pressure fluctuations 
spectra at x/c=0.73 are presented in Fig. 14b. The 
amplitudes of the fundamental are similar, but the 
comparison is not completely unbiased since the mean 
flows are different (Fig. 13). The F^O.5 flow is still 
separated at x/c=0.73 while the F*=l excited flow 
undergoes a pressure recovery. Taking this into 
account, it can be observed (Fig. 14b) that excited flow 
at F*=0.5 is dominated by the excitation frequency, 
while the flow excited at F*=l shows stronger activity 
at higher harmonics. Both frequencies eliminated the 
F + =0.3 mode that was active in the 3D baseline flow. 
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Fig. 14c presents the baseline and excited flow at 
x/c=0.85, where both excited flows undergo a pressure 
recovery. It is interesting to note that 2F for F*=0.5 is 
active, along with F, while for F^l, 3F is active. 
Another feature to note is that F* = 1 decays faster than 
F*=0.5. This trend also continues at the reattached flow 
(i.e. x/c=l.l 1, Fig. 14d and x/c= 1.42, Fig. 14e). The 2F 
of F*^^ and 3F of F f =l continue to play a role in the 
recovery process of the boundary layer downstream of 
reattachment, with the decay rate of 3F (i.e. F*=3) 
being the smallest. Nevertheless, the coherent pressure 
fluctuations of the F*=l excited flow are smaller than 
those of the F*=0.5 excited flow. 

The effect of increasing the excitation level of 
F*=0.5 from <c p >=0.04% to <c p >= 0.3% on the mean 
and fluctuating wall pressures was also presented in 
Fig. 13. The total C p ’ at x/c=0.67 increased by a very 
small fraction of the order of magnitude increase in 
<c M >, but the effect on the C p is significant. To explore 
this effect, <c M > was calculated (Fig. 15a) for F^O.5 at 
the two excitation levels that were used. The main 
effect of the increased level of the cavity pressure 
fluctuations is to increase the magnitude of all spectral 
peaks by about an order of magnitude, besides the 2 nd 
harmonic that increased almost by two orders of 
magnitude. The relative magnitude of 3F and 5F (with 
respect to the logarithmic decay rate as F increases) is 
larger than 2F and 4F. Fig. 15b presents the wall 
pressure fluctuations spectra at x/c=0.73. The two 
mean flows are different (Fig. 13), so the comparison is 
biased. Still the C p ’ of F and 2F are identical, while the 
high <c M > excitation generated very strong activity at 
3F, 4F, 5F and 6F (solid thin line). The level of the 
random pressure fluctuations was significantly reduced 
in the high <c M > excited flow (i.e. at frequencies that 
are not an integer multiple of the excitation frequency). 

Fig. 16a presents the amplitudes of the 
fundamental excitation frequency for the 2D and 3D 
flows when the flow was excited using F^O.5 and 
<c m >=0.3%. The C p ’f. 3 D amplitude was scaled by 
cos 2 (30°). The excellent agreement between the 2D and 
3D scaled C p ’ f indicate the validity of the conventional 
3D scaling for the reattached flow and supports the 
hypothesis that the wave propagates perpendicular to 
the excitation slot. 

Fig. 16b presents the phase velocities for the 
F^O.5 high <c M > excitation in the 2D and 3D flows. 
Without any scaling, both phase velocities start with 
about 0.3U e (x/c~0.7). The trend of the data is very 
similar, regardless of the sweep angle, from x/c=0.67 
to x/c~ 1. For x/c> 1, the 2D phase speed accelerates to 
above 0.8, while the 3D UyUe accelerates to only 0.6, 
which is too small for large coherent structures in 
attached turbulent boundary layers. If we assume that 
the main direction of 3D flow development is along the 
chord, and scale by U e cos(A) to account for it (Fig. 


16b), the agreement between the 2D and 3D scaled 
phase velocities improves. The agreement is within the 
experimental uncertainty, from the x/c location of the 
baseline separation (x/c=0.65) to x/c= 1.25. 

4.3 The Effect of Steady Suction on the Dynamics of 
Flow Reattachment 

It was shown 1 ' 6 that separation control on the present 
model could be performed with comparable efficiency 
using low c M steady suction or periodic excitation (i.e. 
C p <0.2%). Furthermore, it was shown that low levels 
of suction not only shorten the bubble but also increase 
the level of wall pressure fluctuations under the 
controlled bubble. Fig. 17a presents a comparison of 
mean and fluctuating wall pressures measured at high 
Reynolds number in incompressible 2D flow 
(R c =16xl0 6 , M=0.25, x/c=0.64 slot). It can be seen 
that the application of steady suction with c fl =-0.08% 
shortens the bubble by about 0.2c and moves the Cp' max 
upstream by a similar distance. Note that a large C p ’ 
was generated at x/c=0.67 due to the suction. To 
explore the nature of the unsteady effect due to the 
application of steady suction, the wall pressure spectra 
were computed along the model and were than 
normalized by the spectra of the baseline pressure 
fluctuations in order to enhance the visibility of trends 
that are related to the control input. The scaled spectra 
at x/c=0.67 (i.e., immediately downstream of the slot 
and the baseline separation, Fig. 17b) show two 
amplified bands of frequencies that were generated by 
the steady suction. These frequencies are centered on 
F^l.7 and its 3 rd harmonic F*~5 (3F). To explore the 
possibility that the effectiveness of steady suction at 
these conditions is due to the generation of the 
unsteady excitation and not directly through thinning of 
the boundary layer or increasing the shear upstream of 
separation, excitation with F 1 " =1.6 and <c M >=0.09% is 
compared to the steady suction in Figs. 17. The mean 
controlled pressures (Fig. 17a) are similar, with a small 
exception being the separation at a slightly larger C p of 
the suction controlled flow as compared to periodic 
excitation. For x/c>0.7, both controlled pressure 
distributions are identical. The fluctuating wall 
pressures are also similar, with the flow excited at 
F + =l.6 having C p ’s that are larger by roughly 25% 
when compared to those generated by steady suction. 
Both C p ’ distributions show a secondary C p r max at 
x/c=0.95 that corresponds to reattachment. The spectra 
of the wall pressure fluctuations (Fig. 17b) of the 
F + =l.6 excited flow shows very similar trends to the 
suction excited flow with three dominant peaks, i.e., 
those at F and 3F correspond to the most unstable 
modes of the suction perturbed flow (x/c=0.67). The 
wall C p spectra at x/c=l.ll (i.e. downstream of 
reattachment, Fig. 17b) are very similar for both 
control methods, reducing the level of the random 
pressure fluctuations for F + <1. The flow excited at 
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F*=l.6 shows the highest relative fluctuations level at 
3F, in agreement with the 2D and 3D controlled flow 
(Figs lOe and I4e). The possibility that the spectral 
peaks seen for the suction excited flow are generated 
by a cavity resonance was checked and eliminated. The 
resemblance between the effects of low C M excitation 
and suction indicates again the possibility of enhancing 
the control effectiveness by combining periodic 
excitation and steady suction for separation control 1,10 . 

5. Conclusions 

The dynamics of active control of flow separation on a 
Glauert-Goldschmied type aerodynamic surface were 
investigated experimentally. Most of the results 
presented in this paper are for lower Reynolds numbers 
than those included in two previous publications, but 
the effect of the Reynolds number on the flow was 
shown to be small due to the inviscid nature of the 
curvature induced separation. 

Evidence was found for the existence of a sub- 
harmonic resonance in the baseline 2D separated shear- 
layer. A reduction of 35% in the momentum thickness 
of the upstream boundary layer has a negligible effect 
on the mean and fluctuating wall pressures and only a 
weak effect on the pressure spectra. It is speculated 
that flow separation on the present model is caused 
mainly by the local curvature rather than by adverse 
pressure gradient or viscous effects. The level of the 
wall pressure fluctuations in the 3D separated flow is 
significantly higher than in its 2D counterpart. This is a 
manifestation of the additional 3D streamwise vorticity 
that can roll-up to generate discrete unsteady 
streamwise vortices superposed on the spanwise 
vortices. Regardless of the mechanism, the result is 
shortening of the 3D separated flow region. Controlled 
2D data shows that the separated shear layer is most 
receptive to F + =1.5 to 2 and also amplifies these F*’s. 
Non-linear wave resonance between F^O.5 and F^=I 
plays an important part in the reattachment process. 
Downstream of reattachment the low F*’s decay, with 
the decay rate of 3F being the slowest. It was 
demonstrated that the increased efficacy of low F*, 
high <c M > excitation is through generation of higher 
harmonics and non-linear interaction among the several 
excited waves. The swept flow separation is more 
receptive to F*=l excitation than its 2D counterpart, 
and also reduces the random pressure fluctuations more 
effectively. High-level excitation generates stronger 
coherent wave motion at higher F\ and controls the 
random motion more effectively. Conventional swept 
wing scaling works well also for the phase locked 
pressure wave features. 

The effect of weak steady suction is to promote 
wide band unsteady excitation that evolves in a similar 
manner to the evolution of coherent excitation at 
F*=1.6. The most unstable frequency bands for 
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separation control on the present geometry are F*~1.5 
and FV5. 
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Fig. 1 Cross section of the wind tunnel model, with an excitation slot at x/c=0.64, 
as mounted on the tunnel sidewall. 



Fig. 2 Top view of the 2D model. The thick vertical lines indicate the location 

of the end plates and the x symbols indicate the location of the unsteady pressure transducers. 



Fig. 3 Mean and fluctuating wall pressures of the baseline 2D and 3D flow, 
R c = 2.4 and 16xl0 6 , M=0.25, slot x/c=0.64. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

F + 

Fig. 4 Spectra of the wall pressure fluctuations, 

2D baseline flow, R c = 4.2xlO fi , M=0.25, slot x/c=0.64.. 
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5*[mm] 0{mm] H 
Baseline 1.89 1.47 1.29 

mr 1.110.95 1.17 




(a) 

A =30 deg 




Fig. 5 Boundary layer total pressure profiles measured 1.25c 
upstream of the LE with and without boundary layer thinning, 
R, = 4.2xl0 6 , M=0.25. 



-0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 

x/c 

Fig. 6 Mean and fluctuating baseline wall C's, conditions of 




Fig. 7 Spectra of the wall pressure fluctuations measured at x/c=0.73, 2D Fig. 8a-d Spectra of the wall pressure fluctuations for 2D and swept (3D) 

baseline flow, with and without boundary layer thinning, R c = 4.2x1 0 6 , flows measured at x/c=0.73 (a), x/c= 0.85 (b), x/c=l .1 1 (c) and x/c=J .42 

M=0.25. (d). R c = 4.2x1 0 6 , M=0. 25. 
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Fig 9 Mean and fluctuating wail pressures of the baseline and controlled 
2D flow, R c = 4.2x1 0 \ M=0.25, slot x/c=0.64. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
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Fig 10 Spectra of <cji> input (a), and wall pressure fluci 
x/c=0.73 (a), x/c=0.85 (b), x/c=l.l I (c) and x/c=1.42 (d’ 
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Cp'f F + =0.5 


-+-Cq all F + =0.5 
Cp’f F + =0.5 
“®~Cp‘2f F +s =0.5 


-Ud/Ue F + =1 


-Ud/Ue F+-0.5 


U+/Ue F+=0.5 
-a- U<j)/Ue 2F (F + =0.5) 


Fig 11 C p ’ phase locked to the excitation frequency (a), and phase 
velocity (b). Conditions as in Figs. 9 and 10. 


3D Baseline 



-e-3D, F + =0.5 f <c^>=0.04% 
3D, F + =0.5, <Cp>=0.30% 
-A- 3D, F+=1 .0, <Cp>=0.04% 


Fig 12 Total and phase locked C p f at the fundamental excitation 
frequency and its second harmonic (a), and phase velocities (b). 
F + =0.5, <c u > =0.32%, R. = 4.2x1 0 \ M=0.25, slot x/c=0.64. 


- - - F + =0.5 <Cp>=0.04% 
F + =1.0 <Cji>=0.04% 


Fig 13 Mean and fluctuating wall pressures of the baseline and controlled 3D ® y 

flow, R = 4.2x1 0 fi , M=0.25, slot x/c= 0.64. in 1Ca e in egen< 
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Fig 14a Spectra of <c M > input, R c = 4.2x1 0 6 , M=0.25, slot x/c=0.64, 
F + indicated in legend. 
















E 


3D x/c=0.85 
3D, x/c=0.85 F + =0.5 
3D,x/c=0.85 F + =1 .0 


F+=0.5 <c u >=0-3% 


F+=0.5 <c,i>=0.04% 


3D, x/c=0.73 F + =0.5 
3D, x/c=0.73 F + =1.0 


3D x/c=0.73 Baseline 
3D x/c=0.73 F + =0.5 <^>=0,30% 
3D x/c=0.73 F+=0.5 «^>=0.04% 


Fig 15 Spectra of <c M > input (a), and wall pressure fluctuations for 3D 
baseline and controlled flows measured at x/c=0.73 (b). R c = 4.2x1 0 6 , 
M=0.25, slot x/c=0.64, F + =0.5, <c M > indicated in legend. 


3D, x/e=1.11 F + =0.5 
3D, x/e=1.11 F + =1.0 



30, x/c=1.42 

3D, x/c=1 .42 F + =0.5 

3D, x/c=1 .42 F+=1.0 


-2D U^/U 


• 3D U^eCo^O) 


Fig 14 Spectra of the wall pressure fluctuations for 3D baseline and controlled flows 
measured at x/c=0.73 (b)> x/c=0.85 (c), x/c=l .1 1 (d) and x/c=l .42 (e). R c = 4.2xlO fi 1 
M=0.25, slot x/c=0.64, <c u >=0.04%, F + indicated in legend. 


Fig 16 Phase locked C p ’ at the fundamental excitation frequency (a), 
and phase velocities (b) for 2D and 3D flows. F + =0.5, <C M > “0,3%, 
R c = 4.2x10* M=0.25, slot x/c=0.64. 
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Fig 17a Mean and fluctuating wall pressures of the baseline and controlled 
2D flow using suction or periodic excitation, R c = I6xl0 6 , M=0.25, slot 
x/c=0.64. 



Fig 17b Wall pressures spectra of the controlled 2D flow normalized by the 
baseline spectra. Same conditions as the data of Fig. 17a. 
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